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Abstract. In the article (Rothleitner and Francis 2011 Metrologia 48 187-195) the 
correction due to the finite speed of light in absolute gravimeters is analyzed from 
the viewpoint of special relativity. The relativistic concepts eventually lead to the two 
classical approaches to the problem: analysis of the beat frequency, and introduction of 
the retarded times. In the first approach, an additional time delay has to be assumed, 
because the frequency of the beam bounced from the accelerated reflector differs at the 
point of reflection from that at the point of interference. The retarded times formalism 
is equivalent to a single Doppler shift, but results in the same correction as the beat 
frequency approach, even though the latter is explicitly combines two Doppler shifts. 
In our comments we discuss these and other problems we found with the suggested 
treatment of the correction. 



1. Introduction 



Correction due to the finite speed of light is one of the most controversial subjects in 
absolute gravimetry. For a long time the inconsistencies in the results published by 
different authors did not get proper attention. It's therefore very satisfying that shortly 
following our paper on this subject [2], another group of authors has also published 
their view on the problem [Tj. In our paper, we uncovered the reasons of discrepancies 
in the results obtained over the years, and found that some currently used formulas 
incompletely account for the initial velocity of the test mass, leading to the error of 
up to 1 /xGal. The findings of our counterparts, however, are a way more significant. 
According to the authors, modern gravimeters overcorrect the influence of the finite 
speed of light by about 4 yuGal [1J. Discovering a systematic error that big would 
be very important, as some intercomparisons of absolute gravimeters, especially those 
with different instrument types involved, show significant discrepancy in the results [3] . 
However, upon scrutinizing the analysis done in [TJ, we had to admit that the suggested 
reasoning not only implements the same misconceptions that flawed the conclusions of 
other authors, but also introduces some new ones. We now explain our opinion by 
considering some old and new misconceptions in more details. 

2. Relativity is irrelevant for the analysis of corrections in absolute 
gravimeters 

The Taylor expansion of the Lorenz factor 
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allows to estimate the effects of special relativity, if ever applicable, to absolute 
gravimeters. The maximum velocity V of the test mass in modern instruments does 
not exceed 2 m/s, so the component ^ is only about 10~ 17 . This value is not only 
much lower than the practically reachable accuracy of 10 -10 [3], but is also below the 
Heisenberg's uncertainty as applied to ballistic gravimeters [5] . The principles of special 
relativity are therefore redundant for the analysis done in [TJ . The facts the authors use 
to analyze the correction, such as the double Doppler shift (formula (9)), or the time 
delay (formula (35)) are direct results of the classical mechanics. 

3. For the accelerated motion, additional time delay is necessary besides 
the double Doppler shift to correctly account for the finite speed of light 

The double Doppler shift formula 



relates the frequency fo of the laser beam directed towards the approaching with the 
velocity V reflector, to the frequency /q of the reflected beam. Implementation of this 
formula to find the beat frequency (fo — /g) is significantly different for uniform or 
accelerated motion of the reflector. 

If the reflector moves with the constant velocity Vg, the frequencies f and f Q ' do 
not change in time (figjlk), therefore the beat frequency, as the authors of [1] point out 
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Figure 1. Frequency of the reflected laser beam for different modes of the reflector's 
motion: 

a — uniform motion — the reflected frequency does not change in time; 
b — accelerated motion — the reflected frequency is higher at the point 
of reflection due to the finite speed of the light wave propagation. 

in the formula (10), for any moment of time is defined by 
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In case of the accelerated motion, the reflected frequency does change in time. The 
frequency corresponding to the double Doppler shift occurs at the moment of the 
reflection r, while the interference occurs at the later moment t (figjlb), and the 
frequency is different at the moments r and t. Ignoring this fact would imply that the 
increased frequency propagates instantaneously, that is, the speed of light is infinite. 
So, to get the beat frequency for the accelerated motion of the reflector, it's not enough 
just to replace Vo with Vq + got at pL we also need to assign the appropriate timing 
to the resulting process. The author's formula (14) for the frequency at the moment of 
reflection should actually be written as 
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so that the beat frequency (/g( r ) — fo) is obtained as if the interference took place at 
the moment of reflection. This beat frequency, mistakenly used for the analysis in [1], 



corresponds to the following acceleration: 

9(r) = ^ Q(tf(r) - /„)) =90 + ^% + g r), (5) 

where A is the laser wavelength. Because the interference is actually happens at the 
beam splitter, the frequency should be taken at the moment t, not r. Expressing r via 
t with 

b-V - g t 2 /2 

r = t (6) 

and substituting it to the right side of Ml), we find the reflected frequency /q (t) at the 
moment of the interference. The acceleration found similar to ^ will then be 

g(t)=g + 3 f(V + g t). (7) 

The above formula represents the disturbed acceleration that should be used to derive 
the correction for the finite speed of light. 

4. The gravimeter data model is essential for deriving corrections 

The expression ([7]) defines the observable acceleration, disturbed by the effects of the 
finite speed of light. If the measurand is calculated in the gravimeter linearly with 
respect to the distance intervals, then any disturbed acceleration 

g(t)=g + Ag(t) (8) 

translates to the measured acceleration according to the formula [6] 

g= I (g + Ag(t))w(t)dt = g + [ Ag{t)w{t)dt, (9) 
Jo Jo 

where Ag(t) is a disturbance acting during the measurement interval T, w(t) is the 

weighting function of the gravimeter. The correction for the disturbance Ag(t) is found 

as 

Ag — — [ Ag{t)w{t)dt. (10) 



Jo 

The weighting function depends on the gravimeter's data model that includes 

• the location of the time-distance pairs along the trajectory (measurement schema); 

• the formula used to calculate the acceleration using the time-distance pairs (working 
formula) . 

Some weighting functions for different data models are shown on the fig|2j The weighting 
function always obeys the unit square property: 

T 

w(t)dt = l, (11) 
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which ensures that the component go of (|9| is estimated bias-free for any linear 
working formula. However, the diversity of weighting functions creates a strong model 
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Figure 2. Weighting functions of absolute gravimeters 

a - model Si = goTf /2, multiple levels equally spaced in time (EST) 

b - model $ = V T Z + g Q T? /2, EST 

c - model Si = S + V T % + g T? /2, EST 

d - 2-level schema 

e - 3- level schema 

f - 4- level schema 

g - the model Si = S + V Ti + g Q Tf/2 fit to 4 levels 
h - the model Si = So + VoTi + goTf /2, multiple levels equally spaced 
in distance (ESD) 

i - the uniform weighting function once believed to belong to the above 
model (h) [5] 



dependency for the correction (10). That's why any inference of corrections can not be 
agnostic to the gravimeter's data model. This means that the model used to calculate 
go should agree with the model used to derive corrections. Unless the agreement is 
observed, the corrections would not adequately compensate the disturbances for which 
they have been derived. An example of such inter-model disagreement can be found 
in the paper |7J. The gravity value is calculated by the least squares fit of the model 
Si = z + V Ti + g T?/2 using multiple equally spaced in time (EST) levels (weighting 
function on fig|2b). For the gradient correction, the continuous least squares fit of the 
simplified model Si = g^Tf /2 was usecjf] which corresponds to the weighting function 
on fig|2a. The finite speed of light correction was based on the 2-level mode|§] with the 
weighting function shown on the figj2]i. As the result, both corrections introduce a bias 
of several /xGals. 

As the corrections can not be obtained without the data models, so the results 
and approaches of different authors can not be compared without specifying the models 
they used. In some cases, special steps have to be taken to achieve the model agreement 

% formula (7-19) of [7] 
§ formula (7-11) of [7] 



and enable the comparison [9]. In [lj, however, both the corrections are derived , and 
the results are compared with no model involvement. Some model information was 
probably supposed to be incorporated in the new "least squares correction" introduced 
in the paper. As the term suggests, the models not using the least squares fhjjj] are left 
aside. There are, however, bigger problems with this correction, which we discuss in the 
following section. 

5. Trajectory disturbance: a number or a process? and the "least square 
correction" 

The disturbance Ag(t) of the observed free falling mass motion ([8]) is a function, while 
the correction ( JlO) ) is a number. Not necessarily the disturbance has to be expressed in 
terms of acceleration — it can as well be expressed as changing coordinate, or velocity, 
but it's always a process in time. Analyzing the speed of light correction, the authors 
of [1] do not make clear distinction between the disturbance and the correction. In the 
formula (24), for example, the disturbed coordinate is called the "trajectory bias" and 
found as definite integral of the beat frequency over the entire measurement interval, 
i.e. as number. On the other hand, in the formula (25) the authors find the second 
derivative of the bias in T, thus treating it as function. Most avidly the "process 
vs number" confusion is revealed in the formulas (60)-(62), where the left sides are 
functions of time Az(t), but the time argument t is missing on the right side. The 
missing t is not a typo, as the expressions are obtained as definite integrals in t over the 
measurement interval T. 

The "least squares correction" introduced for the bias does not resolve its 
process/number dualism. First, the authors give an example related to the vertical 
gravity gradient, saying that if a model without the gradient is fitted to the data 
containing the gradient, the fitted value of go will refer to a certain point between 
the start and the end of the drop. Then the authors generalize the example saying that 
"this kind of correction can be done for any perturbation. " Such heuristic definition of 
the correction by example and analogy explains neither what is corrected, nor why, nor 
how. 

Based on the efforts described in the section 5 of pQ , we can suggest that actual idea 
behind the "least square correction" was to determine how the disturbed acceleration 
translates into the measurement bias in case the least squares fit is applied to calculate 
the result. Using the weighting functions, this approach can indeed be very efficient, 
but requires expressing a disturbance as polynomial process in time, and then replacing 
its coefficients with corresponding averaging coefficients p]. This technique, however, 
does not create any new kind of correction, but just helps to evaluate how the known 
disturbances translate into corrections. 

|| For example, the 2-level model used in the paper [7] to derive the correction 



6. The retarded times, while equivalent to a single, not double Doppler 
shift, still lead to the same correction 

Let's consider the delays experienced by the photons reflected from the test mass at 
the moments t\ and t 2 (fig. [3]). Because of the finite speed of light, the photons reach 
the beam splitter at the later moments t[ and t 2 . Let's find how the time interval 
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Figure 3. Distortion of the measured time intervals due to the finite speed of light 



At' = t 2 — t[ is related to the time interval At = t 2 — t\. We have 

. (12) 
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t' 2 =t 2 + ^, (13) 



so 
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At' = t' 2 -t[ = t 2 -tx - c At ~~' ( 14 ) 

where d is the vertical separation between the positions of the reflector at the moments 
t\ and t% (fig. [3]). For a small d, the observed velocity V = d/At' of the reflector is 
related to its true velocity V = d/At as 

v- A_^ i _-^_-^_-_ i ^ ssV ,A + n ( 15 ) 

f-r i-r I-? ~ v { 1+ c)- (15) 

The observed velocity looks like it has undergone the Doppler shift, even though no shift 
was considered. While not totally coincidental^] the similarity of the formulas does not 
enable us to compare physical models based on the implemented number of Doppler 
shifts and signal delays, because the ways they implemented may vary significantly 
from one model to another. 

Deriving the correction through the beat frequency is based on the interference of 
the direct and the reflected beams, so two Doppler shifts have to be applied to the signal 



Both the Doppler shift and the signal delay are different manifestations of the finiteness of the speed 
of light. 



frequency, plus a single time delay applied to the reflector's velocity, as discussed in the 
section [3j The resulting disturbed acceleration is given by the formula ([7]). 

The retarded times approach uses only one time delay, but applied to the reflectors's 
coordinate. The equivalence of this method to the first one can be established by 
comparing the disturbed accelerations. Assuming S(t) = z + V t + g t 2 /2, we get 



which is equivalent to ([7]). This proves that the retarded times approach yields the same 
correction as the beat frequency one. 

7. Miscellaneous notes 

Title of the paper. The second-order Doppler shift is a synonym for the relativistic 
Doppler shift |10| . which is irrelevant for absolute gravimeters. In pQ the term 
is actually used for the second order component of the expansion of the regular 
Doppler shift. 

Disturbance magnitude. The formulas (60) - (62) of pQ disagree in their orders of 
magnitude. Only the formula (61) has the magnitude of a disturbance, while two 
other formulas represent entire disturbed trajectory. 

The correction by Kuroda & Mio. The authors' comments on the way the correction 
is obtained in the paper [TT] are not supported by the content of [11]. Nowhere in 
[TT] we find that the reflected beam frequency is obtained using a single Doppler 
shift, as the authors of pQ claim. The factor 2 in the formulas (5) and (5') of [TT] 
is a legitimate consequence of the correct use of the double Doppler shift. 

8. Conclusions 

(i) For the test mass velocities used in absolute gravimeters, the influence of the 
relativistic effects on the measured gravity is below the discernable level and are 
irrelevant for the analysis of the corrections. 

(ii) The frequency of the beam bounced from the accelerated reflector changes in 
time and distance, causing the frequency at the reflector be different from that 
at the beam splitter, so an additional time delay is needed to correctly deduce the 
acceleration disturbance from the beat frequency. 

(iii) The formulas similar to the Doppler shift appear in the retarded time models 
because both are based on the fmiteness of the speed of light. The beat frequency 
and the retarded time models represent different aspects of the light nature and 
can't be validated against each other based on the number of Doppler shifts 
necessary for their implementation. Both models are equivalent in terms of the 
disturbed acceleration, which is defined by the formula (l7j) of this publication. 
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